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ABSTRACT 


This thesis presents an analysis of the throughput characteristics for several 
classes of Aloha packet networks. Specifically, the throughput for variable packet 
length Aloha utilizing multiple power levels to induce receiver capture is derived. 
The results are extended to an analysis of a selective-repeat Automatic Repeat 
Request (ARQ) Aloha network. Analytical results are presented which indicate 
a significant increase in throughput for a variable packet network implementing 
a random two power level capture scheme. Further research into the area of the 
near/far effect on Aloha networks is included. Improvements in throughput for 
mobile radio Aloha networks which are subject to the near/far effect are presented. 
Tactical Command, Control, and Communications (C3) systems of the future will 
rely heavily on Aloha ground mobile data networks. The incorporation of power 
capture and tle near/far effect into future tactical networks will result in improved 
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I. INTRODUCTION 


United States Department of Defense (DOD) Tactical Command, Control, and 
Communications (C3) systems of the future will rely heavily on ground mobile 
data networks. These networks must provide high reliability and connectivity in an 
environment in which enemy counter-C3 assets will attempt to disrupt and deny 
friendly C3 networks. A great deal of interest has been focused on Aloha network 
architectures for future mobile tactical data systems due to their simplicity and 
robustness. 

Research on Aloha has been conducted since the early 1970s [Ref. 1:p. 253]. Most 
of this work. however, has utilized idealized chaunel models. While this has pro- 
duced a wealth of information about the general characteristics of Aloha networks. 
it remains critical that real world conditions be taken into account. Factors which 
affect Aloha channels in the real tactical environment include fading, Near/Far ef- 
fect. power capture (both natural and induced), and jamming. This thesis presents 
analysis of Near/Far effect and power capture on channels which employ the Aloha 
network discipline. 

Chapters II and III] concentrate on the effects of induced power capture when 
emploved as part of a vaniable packet Aloha network. Chapter IV presents an anal- 
ysis of the ramifications of Near/Far effect on three different Aloha configurations. 


This work is geared specifically toward ground mobile data network applications. 


Il. POWER CAPTURE ALOHA 


Aloha packet networks have been utilized for some time in data communications 
systems which are composed of many uncoordinated users, each with low data rate 
requirements. Recent emphasis on military and commercial ground mobile data 
networks will increase the number of systems which employ variations of the Aloha 
protocol as a channel discipline. Although many Aloha systems have been effec- 
tively employed to date, several variations on the basic Aloha scheme stil] remain 
unstudied. The preponderance of work on Aloha network analysis has focused on an 
idealized theoretical radio channel, thereby neglecting the effect of real world chan- 
nel conditions on system performance. It is particularly important to mcorporate 
channel characteristics ito Aloha analysis because, in contradiction to mtuition, 
the overall throughput of an Aloha channel can actually increase under certain con- 
trived or natural channel conditions. A specific example of this is the case of power 
capture Aloha. 

Another area of neglect in terms of Aloha network analysis is that of variable 
packet Aloha. Although model results have shown that fixed packet and slotted 
Aloha both exhibit better throughput characteristics than variable packet. the real- 
ity of bit stuffing logic in fact reduces the actual throughput of these disciplines. In 
contrast, variable packet Aloha can be shown to be a simple and efficient network 
disciphne. 

This chapter introduces and develops the throughput characteristics of a variable 
packet Aloha network which incorporates the power capture phenomenon. For this 
analvsis. only created capture 1s discussed, although the results can easily be adapted 


to accommodate natural capture due to fading. 


Created capture effects occur when groups of users are assigned different transmit- 
ter power levels in an effort to create priority classes, or when all users dynamically 
select a random transmit power in order to increase channel throughput. Previ- 
ous analysis of power capture in slotted and fixed packet networks has shown a 
dramatic increase in channel throughput over the idealized channel models [Ref. 2]. 
This chapter investigates the effect of created power capture with two random power 


levels on an Aloha network using variable length packets. 


A. SYSTEM DYNAMICS 


In this analysis, it is assumed that there is an infinite user population and that 
the channel traffic rate 1s Poisson with parameter g packets per second. Assume 
newly created packets are of variable length x (bits) and are independently and 
identically distributed according to the probability density function a(z) with mean 
7. Within this model, all channel packets are indistinguishable. That is to say that 
newly arrived packets and retransmitted packets have their length x redrawn afresh 
from a(a) for every packet. This is an analytical device known as the independence 
assumption. 

The probability that m packets interfere (overlap) with an arbitrary tagged packet 


during its transmission is given by 


Pr{n} = lole tw)" —o(e+u) (2.1) 


n! 
where z is the tagged packet length, u represents the random length of a preceding 
packet, and n is the number of interfering packets [Ref. 3]. For the case of n = 0, 
this can be interpreted as the probability that the tagged packet encounters no 


overlap from preceding packets and the next interarrival is larger than x. Given the 


case where no interferers can be tolerated by a tagged packet, equation (2.1) can be 
viewed as the probability of successful tagged packet transmission. 

This analysis, however, considers the case when a receiver can be captured by a 
tagged packet in the presence of interfering packets if the tagged packet signal power 
exceeds the joint signal power of all interfering packets by some specific capture 
threshold, +o. If at any time during the tagged packet interval, the ratio of tagged 
packet power to joint interfering power falls below yo, the tagged packet is considered 
to be destroyed. 

Unlike slotted Aloha where the number of interferers is constant over the entire 
tagged packet transmission interval (Figure 2.1(a)), or fixed packet Aloha where the 
number of interferers is a stochastic process made up of early and late interfering 
packets (Figure 2.1(b)). the number of interferers present in a variable packet Aloha 
channel represents a stochastic process which 1s driven by four classes of interferers. 
Figure 2.]1(c) depicts the four possible interferers types. The terms early and late 
refer to the begin and end time relationship between an interferer and the tagged 
packet. For instance, the interferer identified as Early-Early in Figure 2.1(c) cor- 
responds to any interferer which begins before the tagged packet (early) and ends 
before the tagged packet (early). 

The dynamics of the interfering stochastic process is as follows. We note that the 
arrival and departure of n interfering packets partitions the tagged packet interval 
into several non- overlapping intervals of random length. The exact number of these 
intervals, given a known number of interferers and their corresponding interferer 


type, is found by the following relationship. 


Total Intervals = (number of Late-Late packets) 


+ (number of Early-Early packets) 


7 raterLate 
= Siot —_ 


Late- 
Ear] 
Early-Late 
oNONUEET SMe Fixed Packet Variable Packet 
Interferers Interferers interferers 
(a) (b) (c) 


Figure 2.1: Aloha Interfering Packets 


+ 2 x (number of Late-Early packets) + 1 


Relating this equation to Figure 2.1(c), it is evident, for instance, that a single Early- 
Early or Late-Late packet would partition the tagged packet into 1+ 1 intervals. One 
interval occurs during the absence of the interferer; the other during the presence 
of an interferer. In a different case, a single Late- Early would partition the tagged 
packet into 2+ | intervals. ‘Two of the partitions would have no interferers present: 
the other partition would occur during the period of the interfering packet. ‘The 
Early-Late packet covers the entire duration of the tagged packet and therefore 
does not provide any partitioning. 

Each interval has associated with it a total number of interferers /(t;), where ¢; 1s 
the start time of the a’th partition during the tagged packet interval. Hach unique 
sequence of total interferers per partition (/(t;) for all 7) shall be called a realization. 


Figure 2.2 depicts a specific realization of interfering packets for a case where n = 3. 





* of Inteferers 


Figure 2.2: An Example of an Interference Realization 


For this specific realization, the tagged packet total intervals is 5. The total number 
(/(t,)) of interferers present during each interval is indicated at the bottom of the 


figure. 


B. MAXIMUM NUMBER OF INTERFERERS 

A computer program (Appendix A) was written to generate a list of every possible 
interfering realization for a given 2. From this list of all possible realizations, the 
number of realizations (C’;) in which the maximum number of interferers present 
at. any time during the tagged packet interval does not exceed sonie number 7 was 
tallied. ‘Table 2.1 gives the results for n = 1,2,3,4,5. 

An example of how one would generate a C;(n) value is as follows. Consider 
the simple case of n = 1. In this case, only four realizations are possible, nainely 
the interfering packet can be an Early-Early, Early-Late, Late-Early, or Late- Late. 


Continuing with this simple example, it can be seen that given there is 1 known 


0.0208 
0.0061 
0.0018 





TABLE 2.1: Maximum Number of Interfering Packets Given n 


interferer, there is no realization where the maximum number of simultaneous in- 
terferers during the entire duration of the realization is less than or equal to 7 = 0. 
In fact, in this example, the maximum number of simultaneous interferers for each 
of the four possible realizations is 1. Therefore, the value of C,(n) for n = 1 1s 4 
as seen in Table 2.1. For the case of n = 2, there are 14 combinations of 4 packet 
types taken 2 at a time. These 14 pairs of interfering packets generate 14 unique 
realizations. Of these 14 realizatious, 4 have a maxunum simultaneous number of 
interferers of 1, and 10 have a maximum of 2 interferers. The program in Appendix 
A is used to automate the process of generating combinations of packets, coustruct- 
ing the resulting realizations, and determining the maximuin nuinber of interferers 
per realization. ‘Table 2.1 is the result of the program. 

Each C;(n) gives the number of realizations in which the maximum number of 


interferers at any time equals 7, given the total number of interfering packets is 


=~] 


known to ben. The function /7 is simply the total number of all possible realizations 
given n packets are known to interfere. The function f(m) is the inverse of J7 and 


is simply used as a shorthand notation for 1/J7. 


C. VARIABLE PACKET ALOHA WITH RANDOM POWER LEV- 
ELS 


Data packets arrive at the receiver after being spatially attenuated. Each packet 
is transmitted with one of two normalized power levels given by the set Q = {1, M}. 
Here, Af is some integer multiple of the threshold yo times the lower power level 
P = 1 (normalized). that is Af = N+ where N is an integer greater than 1. Each 
user selects a transmit power level for each individual packet from the set 92 with 
equal probability ; in order to avoid any class prioritization among user groups. 


The higher power level A/ is chosen according to the relation 


(N+ 1) > M > Noo (2a 


where iV > 1 is an integer and 79 is the power capture threshold of the receiver. 
The tagged packet, arriving at the receiver with power P; € {1, Af} may capture 


the receiver given a realization of n interfering packets if and only if 


I(t,) 
P, > yomaz[)_ P;] (225) 


j=l 
for all intervals (2’s) in the realization where P, is the power level of packet 7. 


For a tagged packet to capture the receiver, two events must occur. First, the 
tagged packet must be of power Af while none of the interferers are of power M. 


Call this event A. then 





IAS Ae = = 


(n+1) 


— 


[ 


o 6) 


The second condition for capture is that the tagged packet must satisfy (2.3) where 


P, = M and all P; = 1. Given a specific realization, this equates to 


N > maz(I(t;)) (Zea) 


The value N is therefore the maximum number of interferers that can be tolerated 
at any instant of time during the tagged packets transmission interval. If event B 


denotes the tagged packet power satisfying (2.5), then 


- 
eA ola t(n) ). Ona) (2.6) 
j=0 
where C;(7) is the number of realizations, given n interfering packets, in which the 
maximum number of simultaneous interferers equals 7. Table 2.1 provides a list of 
C(n)’s and the inverse f(r) of the total number of possible realizations given n 
interfering packets. 
By using equations (2.4) and (2.6) the conditional capture probability Pr{capture|n} 


is obtained. 


Prictpiureian =O = | (an) 
PrVoann cia = be" Pri Aln} Pri Bin} (2s) 
¢ 
= eee) Se C7 (22) (2.9) 
j=0 


The overall probability of receiver capture by a tagged packet is found by mul- 
tiplying the conditional probability in equation (2.9) by the probability that there 
are n packets in the channel, and summing over all n. Equation (2.1) provides the 
probability of n packets, therefore, the probability of receiver capture by a tagged 


packet which shall be called the probability of success is 


lee 


g(x tu 


n! 


N 
Pr{success} = e 9\7t")[] + » Ota (nr) 9(@ + WIP 2G (n) (2.1) 


The conditional channel departure rate in the interval z + u is given by Reference 3 


to be 


o(r.u) = g(x + u)Pr{success} (221g) 


The channel throughput S, in packets per mean packet time, ts related to the channel 


departure rate by 


ane 
= ca if 5¢ (a, uja(x)a(u)drdu (2p) 
Define G to be the attempted packet rate in packets per mean packet time (i.e. 


G = g®). Appendix B contains the detailed solution of the above equation. The 


resulting channel throughput is found to be the following: 


ee : n43) 5 | 
= Eeaye +e f(n)G?*"(G 41)7! rine 2 Cal) (22s) 


Figure 2.3 shows the throughput of variable packet Aloha with random two power 
signal for values of NV = 0,1,2,3. Packet lengths are exponentially distributed and, for 
computational purposes, n 1s taken to be between 1 and 5. The curve corresponding 


to N = 0 is for conventional variable packet Aloha. 
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Figure 2.3: Throughput of Variable Packet Aloha With Power Capture 


Ill. SELECTIVE-REPEAT ALOHA WITH 
CAPTURE 


To continue the study of variable packet Aloha, the results of the previous chapter 
are now applied to a more advanced Aloha network which incorporates selective- 
repeat as a data link control. 

The specific architecture discussed here is a network where remote terminals com- 
municate with a hub station via an Aloha channel in the query-response mode. The 
feedback hub-to-remote link is an asychronous time-division-multiplexed (ATDM) 
channel. Each remote terminal receives the same data stream and selects applicable 
packets by searching for its own address. 

Assume that packet length of newly arrived and retransmitted packets in the 
Aloha channel are independent and identically distributed with probability density 
function a(z) and with mean ¥. In this model, packets are indistinguishable. Both 
newly arrived packets and retransmissions can be pictured as having their length 
redrawn afresh from a(2) prior to each transmission. 

Packets in the ATDM channel are independent and identically distributed with 
density function b(y) and mean 7. 

Variable packets on the Aloha channel are broken up into z numbered and con- 
stant length minipackets prior to transmission. When packet damage occurs, only 
the specific minipackets which are unreadable are retransmitted. Let D(z),z =1,2... 
be the distribution of the minipackets. On the average, each large packet is seg- 


mented into Z minipackets. 


Let the channel attempted packet rate be g packets per second and Zg the channel 
attempted rate in minipackets per second. With m being the fixed length of a 
minipacket, the attempted channel rate in minipackets per minipacket time is mzqg. 

The conditional probability P, that a minipacket encountered no collisions is the 
probability of no overlap from preceding packets and that the next interarrival is 


larger than the minipacket size m. Thus 


P,(z) = e9 tt) (3.1) 


In a capture environment, however, P; must take into account the times when the 
tagged packet captures the receiver in spite of the presence of n interferers. There- 
fore, using equation (2.1) from the previous chapter, the probability of minipacket 


SUCCESS 1S 


oC 


n N 
—g(mz+m ao-—(n g Te an mm Q- 
Paz) = ML + DT F(a > CXn) (3.2) 
: j= 


| 


The couditional probability of successful minipacket transmission 1s 


q(2,y) = Ps(z)[1 — pi(m)][1 — pa(y)] (3.3) 
where p,;(m) is the minipacket error probability due to errors in the Aloha channel 
and p2(y) is the error rate in the ATDM channel. The conditional minipacket 


departure rate in the interval of = + 1 minipackets 1s 


Onl sald! eat Glee) (3.4) 


The average minipacket throughput is given by the following equation: 


Assume that each packet consists of a geometrically distributed number of fixed 
length minipackets, then the distribution D(z), which is the probability that a packet 


consists of z minipackets, is defined as follows. 


D(z) =(1—y4)?"p (3.6) 
where y is the probability that a minipacket is generated. Because 7 = WEG D@ 
can be written as 
1 ee le 
De) a a (3.7) 


Appendix C provides the detailed solution for the throughput equation. The 


resulting rather lengthy set of equations is as follows: 


OO 


o(1 —p,(m)) Ss h(n)w*[F(n)R(n) — BI (3.8) 


7—1 
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oGv*(a — v)(1 — p(m)) 


ne 2(Z + 1)w*(1 — v)? oe 
Pin) = (oe ie (3.10) 
R(n) = Sov" S7(-1) E(k, r) (3.11) 
r—1 k=0 
(r — hk)! a 
Eka) klik +1) (31) 
G = mzZg (3. ly 
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N 
ee 2 eet 0) (3.14) 


va (1 a)e# (3.15) 
w=(1- =) Ce Me) 
k=z+1] (salen) 
Sar | [sale } 


Figure 3.1 shows the throughput for the variable packet Aloha with selective- 
repeat data link control. The parameters are: p,(m) = 0, po(m) = 0. x and y are 


both exponentially distributed, F = ¥, and = = 8 minipackets. 
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Figure 3.1: Throughput of Selective-Repeat Aloha with Power Capture 
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IV. NEAR/FAR EFFECT ON ALOHA 
MODELS 


This final analysis chapter deals with a network phenomenon which is specific to 
mobile packet radio networks that do not incorporate adaptive control of terminal 
transmit power. The near/far effect is a condition which occurs when packets arrive 
at a receiver with different mean power due to differences in transmission distances. 
Under these conditions. a tagged packet may be able to capture the receiver, in spite 
of interfering packets. if the tagged packet power is sufficiently greater than the joint 
power of the interfering packets. Because each packet is spatially attenuated due 
to the near/far effect. all packets have a finite probability of capturing the receiver. 
The throughput characteristics of various Aloha network configurations which are 


subject to the near/far effect are derived and presented. 


A. CELL MODEL 
The mean power w of a packet at a distance r from the transmitter is of the 


general forin 


Se Tee (4.1) 


In the event of ground-wave propagation without shadowing. 


C= P7,Gr7:Grie,H;} 


Pa, 
des 
to 

eS 


where Pr;, Gr;, and Hy; are the transmit power, antenna gain, and antenna height 
(above ground), respectively, of the transmitting mobile terminal. The values Gp 


and Hp are the gain and height above ground of the base station antenna. The 


i 


exponent @ gives the attenuation law for the channel considered (2 < a < 5). In 
the event of UHF propagation, a typical value is a = 4 [Ref. 3:p. 264]. 

Assuming identical mobile terminals and omnidirectional antennas c is normal- 
ized to unity for all users without loss of generality since receiver capture power is 
determined by the ratio of signal powers [Ref. 3:p. 264]. Additionally, the value of 
a in this analysis is taken as 4. 

In this model, signal power depends on the spatial distribution of users in the 
network. Users create a traffic density g(r) within a circular cell around a base 
station receiving node. The density g(r) can be thought of as the normalized packet 
traffic per unit area at a distance r from the base receiver. The packet generation 
rate in a given area depends on the distance to the base station r and is independent 
of direction. Given this cell model, the total traffic rate offered in the network can 


be described as 


Ge 7 a nomad (4.3) 
0 


Figure 4.1 depicts the general cell model used in this analysis. 
The spatial distribution of packet generation is found to be the probability a 


packet is generated within a distance R. That is to say 


Th ah 
eC Re eps ae ie) = | rg(x)dx (4.4) 


Differentiating, the probability density function for packet generation is found to 


be 


fr(r) = arg(r) (4.5) 
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Figure 4.1: Cell Model for Near/Far Effect Analysis 


in Figure 4.2. As seen in Iigure 4.2, g(r) is normalized to G/ar 


to define a probability density fuuction for the packet power w. 


he 


As a means of simplification, users are assuined to be uniformly distributed in a 
circle of radius r,,,, With the base receiving station at the center of the cell as shown 
loi 7:45) 


From this, the equations for the total network traffic and packet generation reduce 


A standard change of variable operation is performed on equations (4.1) and (4.6) 


max 


Figure 4.2: Distribution of Users Within the Cell 


2 2 
f(w) = we (4.8) 


for. 7. Saco: 
At this point, the network cell model has been described with the distribution 


of packet generation given by equation (4.6) and the packet power density function 


given by equation (4.7). The near/far effect is now incorporated into the model. 


B. NEAR/FAR EFFECT 


In a phenomenon which is similar in many respects to the power capture case 
described earlier, the near/far capture occurs when the ratio of a tagged packets 
power to the joint interfering power exceeds some receiver capture threshold. Again, 
the number of interferers during the tagged packets transmission is a stochastic 


process. The capture ratio is described as 


20 


i (4.9) 


where X is the packet power of the tagged packet and Y is the joint interfering 
power of some number n interfering packets. As in previous discussions, the receiver 
is captured if the ratio described by Z exceeds the threshold given by the parameter 
‘Yo: 

The distribution of X is simply the single packet power density function given 
in (4.7). The joint interfering power density Y is dependent on the number of 


interferers and is given as 


fy(yln) = [fr(w)]"* (4.10) 
where n@ denotes the n-fold convolution. For the cases of n = 0 and n = 1, no 
convolutions are required to determine fy(y). For higher numbers of interferers, a 
computer based convolution package is used to generate specific values for fy(y|n). 
Appendix D provides details on the analytical result for the case of n = 1 and the 
numerical results for m= 2.3. 

Utilizing the density functions for X and Y, the distribution of Z can be found 


using the following formula. 


fz(2|n) = ot yt x (yz) fy (y|n)dy fae) 


The capture probability, that is the probability that Z exceeds the receiver thresh- 
old yo given n interferers, is found by integrating fz(z|n) from yo to oo. 

Table 4.1 presents the probability of capture given n interferers for the range 
0O<n <4. Although the case of nm = 1 is in a closed form with parameter 7, all 


subsequent capture probabilities assume a value of 49 = 3. 
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Pr{capture|n} 












20 
0.057 
0.022 
0.010 


—2/ev 
1, / 







TABLE 4.1: Capture Probabilities 


The limited list of capture probabilities listed in ‘Table 4.1 is sufficient to provide 
excellent resolution of the throughput characteristics for the Aloha models discussed 
in the following sections. The probability of capture for cases where there are greater 


than 4 interfering packets 1s so small as to be considered insignificant. 


C. ALOHA NETWORK MODELS 
Three common Alolia models are analyzed with respect to the improvement in 
overall network throughput gamed from the near/far effect. In each case, the general 


definition of throughput is 


© = \Gigagetccess} (a) 


which in a capture environment can be viewed as 


INO 
iO 


S=G)> Pr{capture|n}Pr{n} (4.13) 


n=O 


The following sections use this approach to determine the network throughput for 
slotted and fixed packet Aloha networks. Variable packet Aloha is analyzed in the 
same general manner as was done in Chapter I for power capture. 
1. Slotted Aloha 
For slotted Aloha. the probability that n packets are generated during a 


given packet time is given by the Poisson Distribution [Ref. l:p. 255]: 


Gre-G 
i! 





Br, = (4.14) 


Therefore. the general form for the throughput in a slotted Aloha network with the 


near/far effect is 


Gr eG 


n! 


SG > Peicapiureln| 


n= 





(Ea) 
Using the Pr{capture|n} values in Table 4.1. the throughput S$ is calculated and is 
shown 1n Figure 4.3. 
2. Fixed Packet Aloha 
The number of packets generated during a given packet time in a fixed 


packet network is given by [Ref. 4:p. 11]: 


(2G)" 96 








OCS ee amar eal: (4.16) 
n! 
resulting in a throughput of 
= 2Gie 
Gaye i ean L e 2G (4.17) 
ae, a 
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Figure 4.3: Slotted Aloha With Near/Far Effect 
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FIXED PACKET ALOHA WITH NEAR/FAR EFFECT 
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Figure 4.4: Fixed Packet Aloha With Near/Far Effect 


Again, using the numerical results for the probability of capture given in ‘Table 
4.1, the throughput is calculated and is depicted in Figure 4.4. 
3. Variable Packet Aloha 
This analysis is very similar to the work presented in Chapter Il} and 
therefore the variable names used here are the sane. First, the probability of success 
for a tagged packet is defined as 


P(z,u) = > Pr{capture|n}e It") [g(x + u)] 


n=O 


(4.18) 


n! 


For this analysis, perfect channel conditions are assumed in order to simplify the 
calculations. The conditional probability of a successful transmission in this case is 


seen to be 


q(x, u) = P(r, u) (4.19) 


The conditional channel departure rate in the interval xz + w is 


O(x,u) = g(x + u)q(z + u) (4.20) 


From this, the channel throughput can be calculated from the following equation: 


il O° Oo 
> =| : o(z,u)a(x)a(u)drdu (4.21) 
2 0d 
Take the case where a(2) and a(u) are exponentially distributed, the above equa- 
tion simplifies to 
: G 


. one) n(n + 2) 5 n+ : =(n-43) . 
S= (Gap +)>_ aa! r{capture|n}G"T "(G+ 1) (4.22) 


Using the numerical results for Pr{capture|n} given in Table 2.1, the throughput 


aoe 


Lj oom | 


for variable packet Aloha with Near/Far effect is shown in Figure 4.5. 


VARTABLE PACKET ALOHA WITH NEAR/FAR EFFECT 


8.18 


0.16 


w 
1 
co 
=) 
oe 
= 
cc 
= 
& 
Ge 
ae 
— 


Ee Ec ae te 
CHANNEL TRAFFIC RATE - G 





Figure 4.5: Throughput for Variable Packet Aloha With Near/Far Effect 
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V. CONCLUSIONS 


This thesis has demonstrated the effects of induced two level power capture on 
variable packet Aloha networks and the effects of the near/far phenomenon on var- 
ious Aloha configurations. In all cases, significant improvements in overall network 
throughput were achieved. 

The development of highly reliable Tactical data systems of the future depend 
on Aloha as a multiaccess scheme for ground mobile terminals. It is critical that 
these networks be allowed to achieve their maximum potential. Given the increased 
reliance on tactical data communications, Aloha networks which provide the maxi- 
mum user throughput possible, even in a jamming environment, must be developed. 

The increases in overall network throughput presented in this thesis indicated 
the improvements in performance which can be realized when power capture and 
the near/far effect are incorporated. As seen in Chapter 2, induced power capture 
increased the network throughput of a variable packet Aloha network from 0.148 
packets per mean packet time in the idealized case to 0.184 in a capture environment. 
Utilizing the selective-repeat data link control with power capture variable packet 
Aloha increased throughput from a maximum of 0.256 to about 0.354. In the final 
analysis chapter, improvements in throughput for slotted, fixed packet and variable 
packet Aloha which incorporated the near/far effect were presented. 

Continued work in the area of Aloha network analysis will be needed in order 
to develop the high capacity, high reliability tactical data networks required in the 
future. Incorporating power capture and the near/far effect into future designs will 


result in better systems analysis and design for these tactical mobile networks. 
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PROGRAM CJ 

Joe McCartin, January 1989 
muoo Pascal — IBM PC 
Bescription: 


This program is designed to generate all 
realizations of interferers given n packets are 
known to interfere. Using the realizations, a 
table of Cj(n) values is compiled. This program 
1s geared specifically toward the purposes of this 
thesis and 1s not a general purpose program. 


2K RK KK KK OK KK KK KK Kk KKK KK KK Kk kK aK RK ok aK ok ok oo ok ok ok ok ok oie oe ok 


program cj (filevar, output); 


const 


Max_intervals = 11; 
matrix_length = 1024; 
@jemacrix_dim = 10; 


type 


realization_vector = array [1..max_intervalis] of 

integer ; 

realization_matrix = array [1..matrix_length] of 
realization_vector; 

eqmvector = array [1..cj}_matrix_dim] of integer; 


Cj_array = array [1..cj_matrix_dim] of cj_vector; 


index = array [1..4] of integer; 
comb_matrix = array [1..220] of index; 


realization_x, 
realization_y:realization_matrix; 


vector:comb_matrix; 
vector_length: integer; 


active_matrix, 
active_matrix_length: integer ; 


eiaematrixve) annays: 
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} 
} 
} 
} 
} 
i 
Ir 
} 
} 
} 
Jr 
ii 
ii 
} 
i 
} 
ii 


init_interferers, 
max_interferers, 
min_interferers, 
num_intervals, 
final_interferers:integer; 


n, 
1, 
j; 
t, 
h, 
runs: integer; 


file_name:string[10]; 
filevar:text; 
outfilevar:text; 


ans:char; 
archive: boolean; 


ORO I I I II I a a ak kK aK aK a ak ak ak 
{ i 
{ INITIALIZE_CJ_MATRIX } 
{ } 
{ Tabulated values of Cj’s are kept in a square matrix } 
{ defined by the dimension cj_matrix_dim. This proc } 
{ initialize this matrix which Js calved Cio. = } 
{ } 
{ 


FOO OC I IGRI I a a ak a ak ak ak ak ak ak ak ak ak ak 


procedure INITIALIZE_CJ_MATRIX; 
var 
a 
j: integer; 


begin 
for 1:11) to ¢)_matnixeaineac 
begin 
for j:> 1 to ¢€jmatrixedimadso 
begin 
Ejomatrix [jd Oe 
end; 
end ; 
end; 
{I OK III I I I I a ak ak ak a a ak ak} 
} 


{ INIT_REALIZATION_MATRIX I 
{ } 
} 


{ In this program, realizations are generated one 
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meaeaterval at a time until the entire realization is } 
{ constructed. For each interval, a matrix of all :; 
{ known realizations is read. Additions to these } 
{ partial realizations are then placed in a different } 
{ matrix. This proc initialize both of these } 

jeematrices. } 
{ } 
{GC OI dd Rk kk kc} 


procedure INIT_REALIZATION_MATRIX; 
var 
1, 
j: integer; 


begin 
writeln (’INITIALIZING REALIZATION MATRICES ”); 
writeln; 
for 1 := 1 to matrix_length do 
begin 
for }) := 1 to max_intervals do 
begin 
realization_xli,j] 
realization_y[i,j] 
end; 


| 
OO 


end; 


merceln (7 INITIALIZAT@ON COMPLETE’) - 
writeln; 
end; 


{OOO IO I GK a ak a ak a aka ak} 
it } 
me COPY _X_TO_Y } 
{ } 
{ This procedure copies the contents of realization , 
weematrax Xk into realization matrix Y. ; 

} 

} 


{ 


{OO I I I I I a a ok kk 2 


procedure COPY_X_TO_Y ( x ,y:integer) ; 
var 
1:integer; 
begin 
for 1 := 1 to max_intervals do 


realization_yly,i] := realization_x[x,i]; 
end; 


{IO I a a a a a aK 


a 


COPY2y_ 10g» 


matrix Y¥ into realization Matrix. 


{ y 
{ } 
{ i 
{ This procedure copies the contents of realization }: 
{ } 
{ } 
GOO OO OOOO OOO OO ORR a a aR aK aK aK aK a } 
procedure COPY_Y_TO_X (y,x:integer) ; 
var 
1:integer; 


begin 
for 1 := 1 to max_intervals do 
realization_x[x,i] := realization_yly,i]; 
end ; 


OOOO OO OOOO OR I RR KK ak ak ak ak ak ak ak ak ak} 


PRINT_MATRIX 


{ } 
iL 
{ } 
{ Print_matrix prints out the most recently updated } 
{ realization matrix (X or Y). If the flag MATRIX is } 
{ set to 1, then the REALIZATIONS wectrneostese } 
{ updated and is subsequently printed. If the flag  } 
{ is not 1, the Y matrix is printed. If the : 
{ diagnostic flag ARCHIVE is true, output goes both  } 
{ to the screen and to the data file OUTFILEVAR. } 
{ Else, output is to screen only. }: 
{ t 
{ } 


2K OK oe ok ok ok 2K 26 ok ok OK ok ok ok ok oi ook ok ok ook ok oi ok KK oi oie ac a oe oi oi ok oo oie oo 2k oo ok ok oc ok ok oc oc oe ok 2 
procedure PRINT_MATRIX (matrix, length, width:integer) ; 


var 
1, 
j: integer; 


begin 
for 1 := 1 to length do 
begin 
writeln; 
if archive then 
Wreiteln (ourtilevan e 
for 3 := 1 to widthado 
begin 
if matrix = i then 
begin 
write (realization_x[i,j]); 
if archive then 
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Write Kouttallevan 
realization_x[i,j]); 
end 
else 
begin 
write (realization_y[i,j]); 
me: archive then 
Wiad ce south levware 
realization_yli,j]); 
end; 
end; 
end; 
writeln; 
if archive then 
writeln (outfilevar) ; 
end; 


OO ORG I I a a a: aK ak ak ak ac ak} 


GENERATE_REALIZATIONS 


{ } 
it } 
{ } 
{ This procedure produces a matrix containing all } 
{ possible realizations given an initial number of i; 
jeeinterferers, a maximum number of interferers, and  } 
{ the total number of intervals. A tree method is i; 
{ employed to build up a realization one interval at } 
{ atime. The first interval contains the initial i 
{ number of interferers. This partial realization } 
{ is placed in a realization matrix. On the next i 
{ pass, the partial realization is read and the next } 
{ interval is added to it. These realizations are is 
{ then placed into another realization matrix. This } 
{ process continues until all intervals have been } 
meecriiled in for all realizations. } 
{ : 
{ i 


EK KK 2 KK ok KK 2k kk ok ok KK kK aK ok 2 ko ok aK a 2k kK ok aK ak 2k 2k ok ok oko ok 2k kK OK aK KK 2 2k ok ok ok 


procedure GENERATE_REALIZATIONS (init_interferers, 
max_interferers, 
num_intervals, 
min_interferers: integer ; 
var active_matrix:integer ; 
var active_matrix_length 
:integer) ; 


i ' 


x liststengan: 
y_list_length, 
interval_count:integer; 


begin 


writeln (’GENERATING REALIZATIONS FOR ’ ,max_interferers,’ 
MAX INTERFERERS AND ’ ,num_intervals,’ INTERVALS’); 


{ writeln (init_interferers,’ ’ ,max_interferers,’ 
7 num_intenvals, ° mineinvertercroyes, 
writeln; 


if archive then 
begin 
writeln (outfilevar) ; 
Writeln (outtillevan- 
>REALIZATIONS 2 2K RR RK?) ; 


end; 
realization x[l, 1] °="inivcinvem erers. 
Xx Listewlenath: 4a 1 
interval_count := 1; 
active_matrix := 1; 
active_matrix_length := 1; 


while interval_count < num_intervals do 


begin 
y_list length := 0; 
{ writeln (within x tomy coo ms 
for 1 := 1 to xilistllengt nade 
begin 


if realization «|i. interval -countile 
max_interferers then 
begin 
y_list_length := y_list_length + 1; 
COpy_Xx_to_y (agyolisre bengal 


realization_ylLy_list_length, interval_count+1] 
>= realization_x[i,interval_count]+1; 
end; 


if realization_x[i,interval_count] > 
min_interferers then 
begin 
ylist length yolistoVerngtnai. 
copy_x_to_y(i,y_list_length) ; 


realization_yly_list_length,interval_count+1] 


-= realization_x[i,interval_count]-1; 
end; 


5-1 


end; 
INbERValeeGoUuNnt >= interval2count + 1; 
aCtIVesMaruix .= 2: 
active_matrix_length := y_list_length; 


x_list_length := 0; 


af interval _count™< niim_ifttervals tien 
begin 
for Ke= ltcomy_] isiee] enpitiwdo 
begin 
“f wma teln s@ewithda, y temx loop’) ;} 
if realization_ylk,interval_count] 
< max_interferers then 
begin 

X_list_length := x_list_length 

sie 
COpyeyotoex (kK jsqaeist length); 


realization_x[x_list_length, interval_count+1] 
:= realization_y[k,interval_count]+1; 
end; 
if realization_y[k,interval_count] 
> mMin_interferers then 
begin 
x_list_length := x_list_length 
aan 
Copy oy voex(k, x list length); 


realization_x([x_list_length, interval_count+1] 
:= yrealization_y[k, interval_count]-1; 


end; 
end; 
interval_count := interval “count + 1; 
active_matrix := 1; 
active_matrix_length := x_list_length; 


end; 


end ; 

writeln (’There Were ’,active_matrix_length,’ 
Reailazations’ ): 

if archive then 


begin 
wWriteln (outfilevar,’There Were 
? ,active_matrix_length,’ Realizations’) ; ena, 


1f active_matrix = 1 then 
print_matrix (1,active_matrix_length, 
invervea lL count) 
else 
print_matrix (2,active_matrix_length, 
interval coum . 


on 


end ; 


OOOO OOO OO I OI I aK a aK I kK ak a a a ak 


PLEGER 


{ } 
{ } 
{ This procedure accepts a realization matrix which } 
{ contains all possible realizations given a specific } 
{ combination of interferers. Each realization in } 
{ the matrix is checked for the invalid condition of } 
{ overcount. An overcount occurs when a realization } 
{ has indicated that a certain number of interferers } 
{ has departed but yet somewhere in the realization  } 
{ the interferer count goes beyond the count of the  } 
{ remaining interferers. Realizations which are } 
{ filtered out are discarded. All good realizations } 
{ are placed into the currently inactive realization } 
{ matrix. J 
* I; 
4. is 


oC tteCtrereeCeCCCCCeCeCCeeeececececeecececerecerce+ err Ls S22 22 25 


procedure FILTER (var length:integer; intervals, 
INCE? remo ee sem. 
var active_matrix: integer) ; 


var 
1, 
j; 
P; 
max, 
num_final_interferers: integer; 


invalid_realization: boolean; 


begin 
writeln; 
writeln (’FILTERING THE REALIZATION MATRIX ’); 
writeln; 


i= ale 


if active_matrix = 1 then 
begin 
active_matrix := 2; 
for 1 := 1 to length do 
begin 
invalid_realization := false; 
num_final_interferers := 
realization_x[i,intervals] ; 
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tienUunetinaimeerterers — anberterers then 


begin 
if intervals > 2 then 
begin 
max := max_interferers; 
fer 9p -—- ero (intervals - 1) do 
begin 
iterealszatienex (i. p} > 
realization_x[i,pti] then 
Max <= Max = 1; 
if realization_x[i,p+1] > max 
then 
invalid_realization 
>= true; 
end; 
end; 
if invalid_realization = false then 
begin 
Gepyax toy (1,3); 
{ print_matrix(2,))3} 
Dee 
end; 
end; 
end; 
end 
else 
begin 
active_matrix := 1; 
for 1 := 1 to length do 
begin 


invalid_realization := false; 

num_final_interferers := 
realization_y[1,intervals] ; 

if num_final_interferers = interferers then 


begin 
if intervals > 2 then 
begin 
max := max_interferers; 
for p 2= 1 fo (intervals — 1) do 
begin 
if realization_yli,p] > 
realization_y{i,pt+i] then 
max := max - 1; 
if realization_y[i,pt+i] > max 
then 
invalid_realization 
:= true; 
end; 
end ; 


if invalid_realization = false then 
begin 


St 


COPy_Y_to=x Gan 
{ print_mat rian, 
]. t= Jae 


end ; 
end; 
writeln; 
writeln (’THERE WERE ’,j}-1,’ VALID REALIZATIONS FOUND’) ; 
writeln; 
if archive then 
begin 
Writelneouttrileyan) = 
writeln (outfilevar,’?FILTER ** >> RK ARK?) | 
writeln (outfilevar,’There Were ’,j-1,’ Valid 
Realizations’ ); 
end ; 
lengths = 7)>i, 
print_matrix (active_matrix,length,intervals) ; 
writeln; 
if archive then 
writeln Wouttiltevar). 
end; 


{OO IO IO IO I I I I I IK IK aK aK a ak ak ak} 
GET Cc 


{ I 
{ } 
“t } 
{ GET_Cj reads in all valid realizations individually } 
{ and determines the maximum number of interferers } 
{ in each realization. A tally is kept of how many  } 
{ realizations have a given maximum number of } 
{ interferers. } 
{ } 
tl } 
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procedure get_cj (active_matrix, active_matrix_length, 
n: integer) ; 


var 
10 
as 
max: integer ; 
begin 
writeln (’FINDING Cis FOR TA REAR ZATION Ve ne 
PACKETS’); 


if archive then 
begin 


Wwrlvetnmcomutfilevar) - 
writeln (outfilevar,’CJ CALCULATION ****#*#«?); 
WEitelnecoutt11eVar) ; 


end; 
mor 1 := 1 to active_matrix_length do 
begin 
if active_matrix = 1 then 
begin 
Nae realization x li. 1); 
for j) := 2 to max_intervals do 
begin 
if max < realization_x[i,j] then 
max := realization_x[i,j]; 
end; 
cj_matrix [n,max] := cj_matrix [n,max] + 1; 
if archive then 
begin 
WRibelnetouttiveyvan) = 
writeln (outfilevar, ’MAX FOR THIS PASS 
1S ? max); 
writeln (outfilevar, ’?CUMULATIVE MAX FOR 
=) lant mee erncibetoc mn, max | ) ; 
end; 
end 
else 
begin 
max := realization_yli,1]; 
for j:= 2 to max_intervals do 
begin 
if max < realization_yli,j] then 
max := realization_y[i,j]; 
end; 
cj_matrix[n,max] := cj_matrix[n,max] + 1; 
if archive then 
begin 
Werte ln oubt thevan). 
writeln (outfilevar, ’MAX FOR THIS PASS 
iSe7 sma) 
writeln (outfilevar, ’CUMULATIVE MAX FOR 
[=">max, ° PKI Ssuan,’ 1S 
we) Wat mieainemasc 
end; 
end; 
end; 
end; 


OOOO OK kak } 


{ } 
{ PRINT_CJ } 


JY 


iL 


i 
{ This procedure prints out the matrix of Cj values } 
{ constructed by the procedure GET_C}. } 
{ I 
jr 
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procedure PRINT_CJ; 


var 
1, 
j: integer; 


begin 
writen (249 2 <3 94 een6 Ses FO oele 
for 12=91 Go Gjlmatrixfaan de 
begin 
Write (i. ee 
if archive then 
Write Coutiilevon sie oe 
for*j:=\1°to ¢cjlmatrixedim co 
begin 
wrate(? 7 yejimatieie lee 
if archive then 
write (outfilevar,’ ’,cj_matrix(i,j]); 
end; 
writeln; 
1f archive then 
Writeln Couttilevan = 
end; 
end; 


ORR I kkk kk a a a k kok kak 
COMB 


i I 
t } 
2 ii 
{ Given a total number of interferers, every I 
{ combination of interferer packet type must be } 
{ determined. This procedure accepts the total } 
{ number of interferers n and generates every i; 
{ combination of the four interferer types which i 
{ add up to the total n. This is basically an tf 
{ enumeration of the number of combinations of 4 ty 
{ things taken n at a time. This algorithm is } 
{ derived from Reference 7. : 
{ } 
{ } 


2K 2 2K 2 OK 2K IK 2K 2k 2K 2K 2K 2K 2K 2K 2K oe ok 2k oo 2K aK KK aK 2K IK 2K 2 2 2K KK 2K KK EK 2K 2K 2k aK OK 2K OK OK 2k ok 


procedure COMB (n,k:integer; var r:index; var mtc:boolean) ; 
label “10, iiibior 20. 30531. 
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val 


1 
Ji, 
W, 
xX, 
y> 
Zz, 
l:integer; 
begin 
HO): vf mtc = true then 
begin 
goto 20; 
end 
else 
begin 
elt \>s= 7; 
t= nN; 
he: = 0: 
el Tome iee:= 2 to k do 
fae — Oe 
15: ieonlk) = n then 
mtc := false 
else 
Nee. =" LGEuUe: 
goto 31; 
ZO: if t > i then 
begin 
hee = 02 
end; 
BO: ho Yeh a 
t e=erih!: 
[Jo || SO aerey 
bi) c=" t = 4s 
r{h+i] := rfh+i] + 1; 
Bove 3), 
end; 
Si: ee |; 
end; 


{OR IOI III II ICI a Ik ka kok a ak ab 


GET_COMBINATIONS 


This procedure controls the execution of the 
procedure COMB and generates the values w, x, y, 
and z. These variables correspond to early-late, 
early-early, late-early, and late-late packets 
respectively. All combinations given a total 


QA A AA AAA 
By YQ YQ YY 


4] 


{ number of interferers n are placed into the matrix } 
{ vector. } 


a } 


{ORI II kk kk kk kk 


procedure GET_COMBINATIONS ( var vector:comb_matrix; 
n:integer; var vector_length: 
integer) ; 


Var 
1, 
j; 
Ww, 
xX, 


ve 
Z: integer; 


mtc:boolean; 


viindex; 
begin 
for 1°5:="1 tov4de 
Viti, = a0. 
1 3;= 07 
mtc := true; 


if archive then 
writeln (outfilevar,’COMBINATIONS 
FO IO): 
while mtc do 
begin 
I? = 0 then 
mtc := false; 
cies ee ce ee 
Como (nl, 4.1 smcele 
for } := 1 to 4 do 
VECLOr (asl = oy Male 
9S ee ll) 
= v[2]; 
sea toll 
:= v4]; 
1! writen (wy, ae ee eee 
if archive then 
begin 
Writeln Coutfilevar,w. 23x, ee eee 
end; 


N< =~ & 
1 


end; 
vector_length := 1; 
if archive then 


writeln 
Coutfilevarr , 722 OOK) ) 5 
end; 
RRO OOOO OO OO a a a a a aK aK aK a aK aK kak} 
MAIN 


{ } 
{ fr 
{ Program execution begins with a query for an If 
{ archive data file. If the archive option is chosen } 
{ all screen I/O will be echoed to the file given. b 
{ The next program query is for the maximum number of } 
{ packets, n. The program will then generate all i; 
{ combinations of interferer packet types for a given } 
{ value of n. Each combination is then fed one by } 
{ one into the GENERATE_REALIZATION procedure where } 
{ all possible realizations are generated. This list } 
{ of realizations is then filtered by FILTER to cull } 
{ out invalid realizations. A matrix of valid | 
{ realizations is fed to GET_Cj in order to generate } 
{ the Cj table. This process is repeated for all } 
{ values of n between 1 and the maximum given by the } 
oeuser. ii 
i fi 
{ } 


3K KK 2K 2K 2K 2K 2K 2K 2K 2K 2K 2K ok OK 2 2 2K 2K 2k 9K 2K 9K OK 2k 9K 2K > 2K 3K 3K 3K 3K > 2K > 3K 3K > 3K OK OK 3K 3K 2K OK OK OK 2k OK ok ok ok 


begin 
archive := false; 


write (’DO YOU WANT A DATA ARCHIVE (y or n) -> ’); 
readin (ans); 
if ans = ’y’ then 
begin 
archive := true; 
write (’ENTER DATA FILE NAME -> ’); 
meadin (file name): 
assign (outfilevar, file_name) ; 
rewErte Coutiilevar) : 
end; 


write (’ENTER THE MAX # OF PKTS PER REALIZATION -> ’); 
meadin (runs) ; 


inie1alize_ejyematrix ; 


fHomem ¢= 1 to runs do 
begin 
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1 to 220 do 
17 Go.4 do 
:= Q; 


age al 
Lona 
vector[i,j] 


if archive then 
begin 
Wiite in 


Coutfilevar , 7 KK 


writeln 
writeln 
writeln 
writeln 
writeln 


writeln 
writeln 
writeln 
writeln 
writeln 
writeln 
writeln 
writeln 


OI AK ak kak kak aka 7) 
writeln (outfilevar,’NUMBER OF PACKETS = ’,n); 
writeln; 

end; 


get_combinations (vector,n,vector_length) ; 
for 1°:= 1 torveetorslenc ride 
begin 
init_interferers 7= veetoul, Mitvect amma 
final_interferers “= vector i, i+ cecomimi 
max_interfeérers <= vector li 9) let vy eetor ine 
vector|1.3] + vectonImed), 
Min interferers == sveencm lea | 
num_intervals := veetorli,2] + 2*vectorli. cles 
Were ene (at .44)] ce ale 


if archive then 
begin 


(outrulevar)- 
(Olitilevan ) = 
Coutfilevar, ?------- 9-9 oer re es 
COutraehe var): 


(outfilevar,’PACKET CONFIGURATION IS ’,vector[i,1i], 
es  VeECtorliy iis: oveeror ies le 
P svercons bat2!]| ie 
(outfilevar) ; 


(outfilevar,’INITIAL PKTS = ’,init_interferers) ; 


(outfilevar,’FINAL PKIS = ?{finaleinterterens 
(outfilevar,’MAX # PKTS = ’ max_interferers)- 
(outfilevar,’MIN # PKTS = ’,min_interferers) ; 
(outfilevar,’NUM INTERVALS= ’,num_intervals) ; 
(outt i levarn) = 
(outfilevar, ?---------- ore oe 

end; 


init_realization_matrix; 


generate_realizations (init_interferers, 
max_interferers, 
num_intervals, min_interferers, 
active_matrix, 
active_matrix_length) ; 
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filter (active_matrix_length, num_intervals, 
final_interferers, 
active_matrix) ; 


get_c] (active_matrix, active matrix length, 
max_interferers) ; 
end; 
end ; 
eat _ Cy ; 
close (outfilevar); 


APPENDIX B 


DERIVATION OF THROUGHPUT EQUATION 
FOR VARIABLE PACKET ALOHA 


This appendix presents the derivation of the throughput equation for variable 
packet Aloha with power capture. 
As seen in Chapter II equation (2.9), the probability of a tagged packet success- 


fully being transmitted 1s 


oc le N 
Dette 9 totes ae f(n)2 (n+) | (2 — 2) X Cnn) (Bale 
rl 
Assign the following identity 
gy" N 
h(n) = f(nj27 ev af 260 (B.2) 


As described in Chapter IIT equations (2.10) and (2.11), the conditional channel 


departure rate and the channel throughput are given by the following: 


O(c sola = ania) (Basi 


2 aa [- o(a, uja(z)a(u)drdu (B.4) 


Utilizing the binomial identity 


Ce ae ( d ea ees (B.5) 


and expanding the exponential term in (B.1), the equation for p,(z, u) can be rewrit- 


ten as 


Pars) e 9 e4 1 7 S h(n) 3 ( uns] (B.6) 
=O 


a=) 


Expansion of ¢(2. uv) in terms of p,(a,u) results in 


Oeum—ag\ ol ule " e-7” 1 a 3 h(n) ‘5 ( a yt (B.7) 


| k=0 


Multiplying produces 


Cem) = |gre *e %" + gue te %| 
Gite eo + gue ae |= ps h(n) SS ( 1 Ayes] (B.8) 
7—! k=0 


Using this definition of @(.7,u) im the equation for S (B.4) produces two separate 
double integrals. The first is equivalent to the equation (27) in Ref. 2. Assuming 
that py(v7) = 0, po(y) = 0 and that the density function a(z) is exponential. this 
first term for S is found in Ref. 2 to be 

€ 


Si =— (era (B.9) 


The remaining double integral remains to be evaluated. This term represents the 


increase in throughput due to power capture. The term is written here as 


— =| / | ¢ nem nen Gem lnem: 7 
2Jo Jo 


. >» h(n) 3 ( , Aus] a(x )a(u)drdu (B.10) 


n= k=0 


Manipulating terms in the above equation produces 
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Further collection of like terms produces 


bie 9 SS h(n) p> ( i {3 aktle-97 a(a)dz ? I ue a(u)du 
ae Ee ( i ; [oem alu)du [ sHe-ra(e)de| (B.12) 
O 0 


Each of the individual integrals above can be solved using the following identity 


(Ref. 5:p. 310]. 


/ rr eda = nly OF) (Bales) 
0 


Additionally, the following substitutions can be made because packet size is as- 
sumed to be exponentially distributed and because the quantity G' is given in pack- 


ets/average packet length. 


a(z) = —e*/* (B.14) 
C= 02 (B.15) 


Combining the integral identity in equation (B.13) with the definition of a(z) 


produces this solution for the first integral term in equation (B.12). 


ean 2) 
/ gk telat es) dy 
O 


_ ED Ley 





ee” 








cin 


—_— 

> 
+ 
= 


The second term in equation (B.12) is evaluated in the same manner as above. 


Oo 
/ u™*e-9a(u)du 
0 


lI 
| 
on 
= 
3 
| 
> 
ao) 
| 
oS 
+ 
ete 
2 
2 
= 


(n — a (Es york +H) (B.17) 
Uu 


Uu 
In this model, all packets come from the same distribution, therefore 7 = u. This 


modifies the above equation. 


54 a bc - (22 = hk)! : =(n—k+1) 2 
i aC tr aie = ==—=0 (G+1) (15) 
0 ie 


The third integral term in equation (B.12) 1s reduced as follows 


a n—-k+1l —gu =. (n a k + Ie) ' —(n—k+2) 
/ ue Maluldu = sea US lym eas (B.19) 
The final integral term is as follows 
a. ki 
/ Beeler )0T rate ijeaet (B.20) 
0 ie 
The total throughput of the system 1s 
So = 5140s 8.2 1) 


Collecting equations (B.9), (B.11), and (B.16)-(B.20) produces 
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= G ee n n! (A ar 1)! : G2) 
le = (G a. 1)7 (n-—k+1) 
=e all k! 
CC a aaa (G + yer (B.22) 


This simplifies to 


G Ga me 
S = (G+15 als 7M (72) > eae 
. an i (n+3) 
+ni(n—k-+ ee a 


Further simplification produces 


a G (Ge ea) ; 
j= (G41) + - eM (7 a — a n+2 (B.24) 
The final summation in the above equation is over the range of k = 0 to co. 


Because there are no / terms within the summation, the entire summation reduces 


to n'n(a + 2). Combining terms produces 
| 


| Ge i | ae (G+ iP eee) . as 
S= (Ga4ip + 59 2 ace SC ylaa a) Se +2)n! (B.25) 


7=0 


a 
—~— 
ed 
ee” 
hh 
” 





Combining terms and utilizing the identity G = g@ produces the following solu- 


tion for the throughput of variable packet Aloha with power capture. 


CO ) N 
: ip SS Al ae ae \GTTN"G a3 1)7 (n+3) S Cu (B.26) 
1 -_ 2=0) 


APPENDIX C 


DERIVATION OF THROUGHPUT EQUATION 
FOR SELECTIVE-REPEAT ALOHA 


This appendix presents the derivation of the throughput equation for the variable 
packet Aloha network with power capture and selective-repeat data link control 
presented in Chapter III. 


From Chapter III equation (3.2), the probability of tagged packet success to be 


(O29 \" N 
= IU Z Thi (ge MZ aE mr oy” 
pa(z) = eramet LES (JOH eB Oslo) @ (C.27) 
il | 


The conditional probability of successful minipacket transmission is given 1n equa- 


fem 3.3) of Chapter II] to be 


q(n,y) = ps(z)[1 — pi(m)]f1 — po(y)] (C.28) 


Also given in Chapter IIJ is the conditional minipacket departure rate in the interval 


of 2 + 1 minipackets. This equation is rewritten here. 


bm(z,y) = mg(z + 1)q(z, y) (G29) 


Equation (3.5) in Chapter III gives the expression for the average minipacket through- 


put and is rewritten here. 


Sm = =—— | bmnlz,y)dy (C.30) 


The equation for minipacket throughput given above is now expanded to produce 


the following: 





Sm = SF ngll = pylons [UL = poly) yey Ue + D(z) 


CO 


mall — pi(myje™® [1 = pay }Olu)dy Yo F(njore» 


Z= i n= 





fez 
z+1 


nei 
Tian Bere 
Agteme + yl" 5° [Csln)] (2 + 1)D(2)e7™9 (C31) 
n! = | 
where D(z) is the probability that a packet consists of z minipackets. 
The first term in the above equation is the same as the equation evaluated in 
Ref. 2 equation (36). The resulting throughput equation for this term is given in 


Ref. 2 equation (37b). Call this part of the total throughput $,. Then from Ref. 2, 


BGO pill) 


‘39 
=(= + 1)w2(1 — v2 ae 
With the following definitions: 
] 
v = we® (C.34) 
G =e nce (C gam 
Fo= / [1 — po(y)]o(y )dy (C.36) 


Let Sp be the remaining term of the total throughput given in equation (C.31). 
By separating terms and noticing that [g(mz+m)]” = (gm)"(z+1)", So is rewritten 


as 








flay GN" Iey(ny}-Se(2 + 1 D(ze-™| (C37) 


n=] ‘ 7=0 z=] 
As stated in Chapter III equation (3.7), assume that D(z) is geometrically dis- 


tributed and is defined as follows: 


D(z) = (1 2 =) (C.38) 


Utilizing the equations for G, D(z), and w the equation for Sz can be written as 











follows: 
G 
i am (1 — mn (m)le7 Za 

- ~cusyilen : . n+l, 2-1 a(S) 

SS f(n)2 ~ S— [C;(n)] (z+1)"" w 2) G3 9)) 

i=) 70 2 ==) 

Assign the following identity: 
j Pec. | a: ye a 
Nt) ae aes ee (C.40) 
. ae 5—O 


At this point. consolidate the exponential terms and divide out two of the w 


terms to equalize the exponents within the summation over z. 


a [1 — py(m)lo SS h(n)w* Be Ju Ly ep? tee) (241) (C.41) 


z£ se ] n=! 2— 
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Assign the following identities: 


a) 
| 
, 
2 
4. 
— 
ee 


(C.42) 


foc= Seely (C.43) 


This produces the following 


OO 


a2 ee 
=n [1 — pi(m) a a (C.44) 


oe 





At this point, modify the summation over k such that the range is from k = 1 
to oo. To do this, subtract off from the total of the summation the value of the 
summation for k = | which is v in this case. Therefore, the above equation is 


changed to 


G 


eS 
Sono iay| 





[1 — p,(m)lo 5 h(n)w~? bs [tv | — : (C.45) 
k=! 


| 


Given that |v] < 1, use the following mathematical identity to solve the above 


equation [Ref. 6: p. 142]. 


ord 7 t r—] ? ae, t 
2. Lty —_ (1 = v) heal + 1)! e8 2al- DEE (C.46) 


Using the above identity, the final form of 52 is given as 





t pal : r — k)t 
{fae (+14 4 1)! ye U- ye | (C4 


The total throughput for selective-repeat Aloha with power capture Is given as 


Ge S (C.48) 


APPENDIX D 


ANALYTICAL AND NUMERICAL RESULTS FOR 
NEAR/FAR EFFECT CAPTURE 


This appendix provides the details of the analytical and numerical methods used 
to determine the values for Pr{capture|n} in Chapter III. First, the derivation of 
Pr{capture|n = 1} is provided. This is followed by a MATHCAD file which was 
used to calculate capture for n = 2 and 3. 

A. ANALYTICAL METHOD 

The case of n = 1 is sufficiently simple to conduct a straightforward mathematical 
derivation of the capture probability given the near/far effect. From Chapter IV the 


following is known: 


L=— (D.49) 
) 
27 
cel ie a (D.50) 
aG 
Fy(yln) = [Sy Pe] (D.51) 
For the case of n = 1, both y and z are defined between the values of 7°, and 


oo. Also from Chapter IV, equation (4.10), the definition of fz(z|n). From this 


CY 


ay Ziven the valid 


equation, it is clear that both yz and y must be greater than r 
ranges for fx(xr) and fy(y|n). This implies that for values of z greater than 1, the 
valid range for y is y > r—°.. Therefore. the equation for fz(z|n) with the proper 


maz’ 


limits of integration 1s 


cn 
Kear t 


fz(=|n = 1) 


Remembering that G = rr? 


fz(z|n = 1) 


mar 8} 


ies ufx(yz)fy(yln = 1)dy 
eZ 

—(2/at1) 27 -Q/ath) gy 
aG 


2 Oo 
(—) y(2/e) 
G 


—a 
Tmax 


&** _-(o/eti) (_ ay" 
aG 4 


-_O 
Tmax 


y~ Alot} dy 





ves us 


~ of 
j 
ae, mar Tmax 


— 1 -eyfo+t) 
CO 


= — 4 _,-@/at) -4/0|* 


(D.52) 


(D.53) 
(D.54) 


(D.55) 


To determine the capture probability from the above distribution. simply inte- 


grate from the capture threshold +9 to infinity. 


Pr{capture|n = 1} 


Prizes s\n =e 


i 1 -@jot1)g- 
yo OQ 


90 


(D.58) 


(D.59) 


(D.60) 


B. NUMERICAL METHODS 


This document is a MATHCAD file used to numerically calculate 
the values of P({capture|n} for values of n from 2 to 3. The 
general approach is to use a numerical convolution to determine 
the n’th convolution of f(y). These convolution results are 
then used in numerical integrations to generate the P(capture|n} 
results. 


Initialize and Define Variables 


Me= 1 ( M corresponds to r ) 
max 
x = 4 
2 
Gin7:M 
pe— 7 ( 2 < D < 32 sets the number of convolution points. ) 


i= 0 ..32°'D- 1 


2;1 +1 

vy, =-— 
1. 22D 

2 
[on 
T x x 
f(t) =2-—t -1¢(10 - t) - fy = e! | 

“:G 

mite) = £(t) 
f(t) corresponds to f(w) in the model 

p(t) = f(t) 


The following equation performs a discrete convolution on the 
functions defined as X and Y. 


1 — a 
conv(X,Y) = -]-sceeel oe), | 
D 


Of 


Step 1. Convolve f(t) with itself to get f(y|n=2) (call this 
C) and plot the resulting functions. 
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Step 2. Convolve f (y|n=2) with f(t) to get f (y|n=3) (called 


D) and plot the relulting function. 
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Step 3. Continue in the above manner for as many cases of 


f(y|n) as needed. 


Cr 
oa) 


Step 4. Calculate f(z|n) using a linear interpolation of 
the numerical value for f(y|n). 


2 2 
-|-+1 30 -|- 
wT x X 
mez — 2: Z y -linterp(v,C,y) dy 
“:G “x 
ANI 
Step 5. Calculate the capture probability given n known 


interferers by integrating f(z|n) (f(z) in this 
example) from the capture threshold (3 in this 
example) to a maximum usefull limit which is 
determined by the range of values for which f(y]|n) 
is valid. In this case, 30 is the upper limit. 


30 
Z i= | f(z) dz 
3 
Z = 0.057 
Step 6. Continue on in the same manner as above to 


calculate the capture probability for different 
values of n. The following equation defines f(z|n) 
for the case where n=3. 


go s= 2° ——- hh Y -linterp(v,D,y) dy 


oe 
ll 


30 
| g(h) dh This equation calculates P(capture | n=3} 
s) 


H = 0.022 


og 
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